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1 Abstract 

The paper discusses a application of multi-dimensional calibration metrics to evaluate pressure 
data from water drop tests of the Max Launch Abort System (MLAS) crew module boilerplate. 
Specifically, three metrics are discussed; 1) a metric to assess the probability of enveloping the 
measured data with the model, 2) a multi-dimensional orthogonality metric to assess model 
adequacy between test and analysis, and 3) a prediction error metric to conduct sensor placement 
to minimize pressure prediction errors. Data from similar (nearly repeated) capsule drop tests 
shows significant variability in the measured pressure responses. When compared to expected 
variability using model predictions, it is demonstrated that the measured variability cannot be 
explained by the model under the current uncertainty assumptions. 

2 Introduction 

During the development of the Apollo command module (CM), many studies were conducted to 
understand the behavior of the CM upon returning to Earth. Because of the limited 
computational capabilities at the time, engineers made extensive use of experimental data to 
supplement the relatively simple analyses. These studies were performed to understand the 
landing behavior under a variety of conditions. Interestingly, Apollo was initially to include a 
land-landing architecture. Thus many of the early studies focused on assessing attenuation 
systems for landing on soil, sand, and clay surfaces [1-3]. However, the land-landing capability 
was abandoned in favor of a water splashdown, similar to Mercury and Gemini [4-5]. Other 
studies in addition to the water landing studies were assessments of the capsule response in the 
event of a pad abort [6]. 
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In 2009 the Orion Project lacked the test data to anchor LS-DYNA™ analytical model 
predictions for the crew module (CM) water landing loads. These loads were the largest 
structural mass design driver and were predicted by modeling the physical interaction of the heat 
shield and side wall structures during water landing. Because these loads drove the design, early 
verification was extremely important. However, water drop tests were not planned until later in 
the program. To reduce the risk associated with using an unverified model for loads, a series of 
tests was undertaken for load verification. Specifically, the NASA Engineering and Safety 
Center (NESC) sponsored an effort to modify an existing full-scale boilerplate CM, referred to as 
the MAX Launch Abort System (MLAS), to use it as a pathfinder for the Orion Project water 
drop tests. Data from these tests were intended to provide an early assessment of the model 
adequacy. 

Recall that Apollo water landing studies built on the extensive experience of the Mercury and 
Gemini programs. A wide range of experimental parameter studies, using subscale test articles, 
was conducted to evaluate the Apollo capsule performance. A primary interest of these studies 
was to understand the effect of landing attitude and longitudinal or horizontal velocity as well as 
heatshield flexibility on capsule response [7-11]. In addition, numerous full-scale tests were 
conducted to assess the capsule performance on more flight-like vehicles, e.g.. References [12- 
14]. Einally, based on knowledge of probable landing conditions, a design envelope was 
established [15]. Although the ORION project has benefitted tremendously from this work, the 
need for verification of computational models has not diminished, and must be undertaken to 
reduce program risk. 

In this paper, data from two capsule water impact tests with the MLAS CM are used to compare 
measured pressures with predictions. However, comparisons are conducted using newly 
developed metrics that consider uncertainty and model adequacy from a multi-dimensional 
viewpoint. In fact the main focus of the paper is to show the applicability of these metrics when 
used with the pressure data. As a byproduct of the application of the metrics, the approach 
fosters the use of sensor placement algorithms to ensure that the data collected are adequate for 
model assessments. Because the placement approach was developed after the tests were 
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completed, only analytieal results are presented for illustration purposes. The paper organization 
is as follows; a test deseription is presented first, followed by ealibration metries, results from the 
applieation of the metries to the test data, analytical sensor plaeement studies, and eoneluding 
remarks. 


3 Test Description 

MLAS test, analysis, and data evaluation were eondueted in two phases over a two-year period. 
The first phase in February 2010 eomprised 18 vertieal drops at four water entry eonditions. 
Tests used a erane at the littoral test faeility at the U. S. Army’s Aberdeen Test Center (ATC) in 
Maryland, to position the eapsule. A picture of the test artiele and a drop test are shown in figure 
(1). On-board instrumentation ineluded accelerometers, inertial measurement units, and heat 
shield pressure sensors and strain gages. In addition, photogrammetric targets on the CM surfaee 
were traeked using high speed eameras to determine the CM position and attitude during the 
water entry. 

For the seeond phase of testing the CM was modified with an expanded number of heat shield 
and side wall pressure sensors and internal strain gages on the primary strueture load path. 
Thirty-eight drops at 21 eonditions were eompleted at ATC in Mareh 2011 with vertieal 
veloeities ranging from 135 to 440 in/see and piteh angles from 0 to 40 degrees. Of all the data 
eollected only tests with impact velocity of 300 in/sec and impact angles 30.6 deg and 21.5 deg 
were selected and used in this study, beeause repeated test data were available. 

The LS-DYNA™ [16] model then generated aeeelerometer, pressure, and strain predictions at 
key test eonditions. The loeations for the aeeelerometers and pressure transducers are shown in 
figure 2. Work here is foeused on pressure data beeause that was the initial foeus of the 
program. 


4 MLAS CM Boilerplate Model Description 

An LS-DYNA finite element struetural model of the MLAS CM strueture, figure (3a), was 
generated by personnel at NASA Ames Researeh Center starting from the CAD geometry. To 
stiffen the test artiele, engineers reinforeed the eore by adding an additional strueture referred to 
as the “egg-erate”, see figure (3b). Most of the strueture (ex. skin, ribs, heat-shield and 
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stiffeners) was modeled using shell elements. However, coneentrated masses were modeled 
using ballast eonnected by rigid link s . Materials models in LS-DYNA used a 
Piecewise Linear Plastieity (material 24) formulation for all structural materials. The heat- 
shield was constructed from ASTM A 516 grade 70 and the crew module was ASTM A572-50, 
except for the top cover assembly, removable hatch cover, CG sensor beam, and bracket which 
were aluminum. The characteristic length of the finite elements making up the model ranged 
from 3” to 0.5”. Water and vacuum volume domains were created and bias-mesh using 8-node 
solid elements with an ALE (Arbitrary-Lagrange-Eulerian) multi-material element formulation 
[16]. Dimensions for the water and vacuum meshed volumes are 240” x 240” x 48” and 240” x 
240” X 12”, respectively, with an element size of 1.0” x 1.0” x 1.0” near the impact region. 
Einally, a penalty coupling method was used for modeling the fluid structure interaction (ESI) 
with a coupling pressure versus penetration curve shown in figure (4). 

5 Pressure Time History Analysis 

The analysis of time history pressure data from water impact tests must be guided by the way the 
data are used during the model calibration. Eor design, peak pressure values are often considered 
as adequate validation metrics used when judging simulation results. Although peak values often 
drives structural design, proper model calibration also needs to capture the spatially dependant 
load variations over the surface of the test article. Because oftentimes high frequency data are 
collected during tests, time independent calibration metrics can provide alternate evaluation tools 
to compare models. These metrics, initially developed for impact problems, are discussed next. 

5.1 Time Domain Calibration Metric-1 

Calibration metrics provide a mathematical construct to assess the fitness of a model in a 
quantitative manner. Work by Oberkampf [17-18] and Schwer [19] set forth scalar statistical 
metrics ideally suited for use with time history data. Metrics in terms of mean, variance, and 
confidence intervals facilitate assessment of experimental data, particularly when probability 
statements are sought. Eor our problem, instead of using response predictions at a particular 
point as a metric, a vector 2-norm (magnitude of vector) of the system response is used instead. 
An important benefit of using this metric is that it provides for a direct measure of multi- 
dimensional closeness of two models. In addition, when tracked as a function of time, closeness 
is quantified at each time step. 
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Because parameter values are uncertain, statistical measures of the metric need to be used to 
conduct assessments. With limited information about parameter variations a uniform distribution 
function was selected and used to represent the parameter uncertainty. This parameter 
uncertainty model is used to create a family of N equally probable parameter vectors with their 
corresponding responses, where N is arbitrarily selected. From the perspective of a user, it is 
important to know the probability of being able to reconcile test with predictions, provided that 
both the model and parameter uncertainty model are appropriate. To this end, let 
Q{t,p) = \\v{t,p)\\^ be a scalar time varying function, in which the vector v, constructed from 
multiple sensors, is used to compute the 2-norm of the response at time t. Furthermore, let 
a{t) = mmQ{t,p)\)Q the minimum value over all parameter variations, and let a{t) = vasLxQ{t,p)hQ 

the maximum value. Using these definitions and N LS-DYNA solutions corresponding to 
equally probable parameter vectors, a calibration metric can be established to bound the 
probability of test values falling outside the analysis bounds as; 

M^=Vxob{git)<QXt)r^QXt)>^)«yJ^ ( 1 . 1 ) 

In here, is the 2-norm of measured experimental responses. Certainly, for cases with data 

from multiple tests, Q^t) mv&i be replaced by its lower QXt) and upper bound 2e(0- The 

arbitrary parameter N controls tightness of the bounds and determines the number of LS-DYNA 
solutions required. 

The use of norms, although convenient, tends to hide the spatial relationships that exist among 
responses at different locations in the model. In order to study the spatial multi-dimensional 
dependency explicitly, a different metric is needed. Before presenting the multi-dimensional 
metric, some needed preliminary time domain decomposition steps are discussed. 

5,2 Orthogonality Metric-2 

Data from either test or analysis is often collected at fixed time intervals for multiple sensors. For 
validation of models to predict loads from pressure pulses, it is preferable to collect data over the 
entire surface being affected. Unfortunately, practical constraints often limit the amount of data 
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collected to smaller regions of the test article. Given n samples of the response sampled at A7 , 
the time history can be decomposed as 
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The rightmost equality of Eq. (1.2) is obtained from singular value decomposition (SVD) of the 
time history data. Data is divided into m measurements stored in Y and q unmeasured values 
stored in 7, . In this decomposition, the basis matrices O e and T e , referred to as 
basis or impact shapes, contains the time independent spatial distribution of sensors, whereas G 
e K'”'" contains time variations, and I contains the non-zeros singular values (cti ■■■ cr„} . 


From Eq. (1.2), one can solve for the unmeasured responses in terms of the measured values as 

= q>(0^0)"‘0^T (1.3) 

After examining Eq. (1.3), it should be apparent that sensor placement plays an important role in 
our ability to predict unmeasured responses Y^. For example, if the matrix O contains 
rows/columns which are linearly dependent (i.e. sensors are placed poorly), Eq. (1.3) becomes 
numerically ill-conditioned and cannot be used to estimate the response . For this reason, 
sensor placement for this class of problem must be guided by the numerical conditioning of Eq. 
(1.3). Furthermore, because some impact shapes contribute more to the overall response than 
others, sensors must be placed to properly recover the most important impact shapes. 
Importance of the z* impact shape is quantified in terms of its contribution to the total response 

as 5,= (1.4) 

/ l=\ 

In this decomposition, the sum total of the impact shape contribution adds up to one. 


With the decomposition approach just discussed, one can now develop a multi-dimensional 
metric to compare models. Specifically, the impact shape vectors computed with Eq. (1.2), 
provide a means to compare models using time invariant metrics like those normally used in 
classical modal test. Impact shapes are ideally suited to compare analytical models with test 
using orthogonality. Numerically, the orthogonality metric is computed as; 

M2=6T (1.5) 
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where O is sized m x h with h measured impaet shapes at m loeations and F sized m x h, are 
shapes eomputed using simulation data. Both O and F have been normalized sueh that $^^0 = i 
and = /. Beeause individual impaet shape veetors are stacked column-wise, the orthogonality 
metric M 2 is sized h x h with diagonal values corresponding to the vector projection numerical 
value. If vectors are identical then their projection equals 1. Consequently, when evaluating 
models, multi-dimensional closeness is judged based on similarity of impact shapes and shape 
contributions. 

5,3 Prediction Error Metric-3 

Two metrics have been presented to assist the model calibration process; Mi quantifies the 
probability of being able to reconcile model with test data, whereas, M 2 addresses model 
closeness from a multi-dimensional viewpoint. Although both M\ and M 2 measure our 
predictive capability between test and analysis, a third metric is needed for this class of problems 
for the sole purpose of conducting sensor placement. This metric considers the error resulting 
from predictions made using data from a reduced set of sensor locations. Intuitively, one would 
expect certain sensor locations to provide a better set of system responses when compared to 
other locations. To quantify this, the prediction error metric is defined as 

M,=(Y^-YyW(Y^-YJ ( 1 . 6 ) 

In here, the quadratic prediction error metric M 3 is weighted by an arbitrary diagonal weighting 
matrix W. This weighting matrix is convenient for cases where the user needs to reject sensor 
locations for other reasons. Since this metric uses only simulation data, it is simple to evaluate 
various candidate sensor sets. To evaluate Eq. (1.6) one uses the exact values for Y^ from 
simulation results and estimated values (using Eq. (1.3)) forF^ . The best sensor set will produce 
the smallest M 3 . In order to obtain a candidate sensor set, combinatorial optimization algorithms 
[20] are ideally suited for the problem. Results later in the paper demonstrate this approach and 
how it compared to the baseline sensor set. 

6 Results using Metrics with Measured Pressure using Metrics 

At this point, our attention is turned to study the pressure data from tests with impact velocity of 
300 in/sec and impact angles 30.6 deg and 21.5 deg. These two cases will be referred to as 
T300-30 and T300-21, respectively. To demonstrate the application of the metrics discussed in 
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Section 5.1, the sensor set was restricted to 33 pressure sensors. These 33 sensors are all within a 
66 inch radius centered at the heat shield impact point, located at x= 779.82, y= -1.25, and z= 
96.91. Although this sensor set is not the full set, many other sensors were left out because they 
did not see pressure values until after 50 msec, which is the time window considered here. For 
the T300-30 condition, data from 7 tests were used, whereas for the T300-21 data from 8 tests 
were used. Recall that metric Mi is used to assess test variation where the maximum and 
minimum singular value bounds are computed using 33 pressure sensors. Figure 5 shows the test 
uncertainty bounds (i.e., Q^(t)andQ^(t) , solid-blue) for T300-30. To construct the tests bounds, 
the magnitude norm using 33 pressure sensors is computed for each of the 7 tests first. Then the 
maximum of all pressure maxima g^(t)and minimum of all pressure minima Q^(t) are collected 

and plotted as a function of time. Note that the observed variation (represented by the difference 
between the upper and lower bounds) from 7 very similar tests is large. Because data from 7 
tests fall within these bounds, it is reasonable to assume that the probability of observing new 
test data outside these bounds is less that 1/8. 

To contrast the uncertainty in the experimental data to predictions from the analytical model, the 
LS-DYNA model was executed multiple times assuming that the impact velocity, pitch, and yaw 
angles were uncertain. Table 1 shows the nominal conditions for T300-30 along with selected 
parameter upper and lower bounds. Uncertainty data reported by the test group was used to 
determine initial parameter bounds but in the end, the bounds used were significantly larger than 
those provided by the test group. With this uncertainty model, parameters are assumed to be 
equally likely between bounds, and as such parameters were varied to generate 28 impact 
conditions near the T300-30 nominal case. As with test data, bounds for predicted pressure data 
are collected from N=28 simulations. Results for the analytical bounds (a lower bound and a 
upper bound) are overlaid to the test data in figure 5; dashed red lines correspond to analysis. 
From metric Mi (shown in figure 5) and the selected uncertainty model, the probability of 
predicting analytical results outside the analytical bounds is less that 1/28. Hence, with the 
current LS-DYNA model and uncertainty models, it is difficult to explain test observations 
outside the analysis bounds. This is in spite of the fact that the parameter bounds used in this 
study were significantly larger than those reported by the test group. Results for metric Mi and 


8 



T300-21 are shown in figure 6 using 8 tests and 35 LS-DYNA runs. Because with this model 
each run takes over 28 hours (using a 4 processor 3.3 GHz desktop computer), the number of LS- 
DYNA runs executed is constrained by available time. Parameters for T300-21 are shown in 
Table 2. As with T300-30 cases, analytical uncertainty bounds do not enveloped the test results 
and therefore, with this model it is difficult to reconcile test with analysis. 

Another way to compare test results to analysis is by using the orthogonality metric M 2 . In here, 
time invariant impact shapes or basis vectors computed using pressure test data are compared to 
analysis, Eq. (1.5). Results for test 1 of the T300-30 series are shown in figure 7a, and results for 
test 7 are shown in 7b. Black squares correspond to impact shapes with high orthogonality 
values between test and analysis. About 9 impact shapes exhibited orthogonality values greater 
than 0.6 for test 7 versus only 5 for test 1. Consequently, results from test 7 (T300-30) are closer 
to the analysis model. Across the top of the orthogonality plot are printed the test impact shape 
contribution to the overall response, whereas the ordinates show values from LS-DYNA. For 
example, from figure 7a, the first test impact shape contributed 14% as compared to 19 % from 
LS-DYNA. Results for the T300-21 are shown in figure 8. Figure 8a compares the 
orthogonality values for test 8 whereas figure 8b compares results for test 4. Clearly, test 4 for 
the T300-21 shows better agreement between test pressures and LS-DYNA. 

7 An Approach for Pressure Sensor Placement 

Because the sensor placement approach in the following was not used during test, it is presented 
last for readers interested in more rigorous methods for sensor placement. When optimal sensor 
placement is performed as part of the pre-test analysis, placement is driven by the test goals. If 
model calibration is the goal, sensor placement must focus on providing information to properly 
evaluate the established calibration metrics. For example, when the orthogonality metric M 2 is 
used, if sensors are not strategically placed, it is difficult to distinguish between basis vectors. 
Moreover, if our goal is to conduct pressure predictions, prediction error must also be considered 
during placement. Our approach uses combinatorial optimization to select a fixed number of 
sensors from a large number of probable locations. To make the problem more tractable, the 
area for sensor placement was confined to a circular sector near the impact point. 
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7.1 Post-Test Pressure Sensor Placement Study 

Placement of pressure sensors for the MLAS CM boiler plate test was eondueted based on best 
engineering practices. Specifieally, sensors were spaced uniformly near the expeeted impaet 
point, with increased spacing as one moved away from impaet. Pressure sensors at 60 loeations, 
as shown in figure 2b, were seleeted and used during test. Beeause of symmetry eonsiderations, 
the bulk of the sensor set loeated during test were eoneentrated on the lower right quadrant of 
figure 2b. From test eonditions, the impact point was expected to be near the z-axis eenter line, 
and it moved along that line depending on the impact angle. For our plaeement study, node 
5017879 with coordinates (779.82, -1.25, 96.91) was seleeted as the impaet center and all 
elements within a radius of 66 inehes were considered for sensor plaeement. This area provided 
for a total of 2148 possible loeations where pressure readings could be reported from 
quadrilateral elements in the region. Figure 9 shows a projection onto the y-z plane of the heat 
shield augmented sensor set (2148 loeation) using blue (+), with red eireles (o) showing the 
pressure sensors loeations used during test. Only 33 of the 60 sensors used during test are within 
the 66 ineh radius used in this study. For comparison purposes, it is assumed that only 33 
sensors can be reloeated to see if the metries ean be improved using an alternate sensor 
plaeement teehnique. 

7.2 Analytical Basis Vectors for MLAS CM Water Impact Simulation 

Because orthogonal deeomposition for this elass of problems is not performed routinely, it is 
important to understand the spatial time independent information eaptured in the basis veetors 
computed using Eq. (1.2). In particular, if characteristic patterns appear from pressure data, 
these patterns can help tremendously in plaeing sensors and more importantly in predieting the 
response at other loeations. For this reason, a set of basis vectors was computed using simulated 
data and all 2148 possible sensor locations. Figure 10 shows contour plots of the heat shield 
region in the y-z plane for the first 4 basis veetors using simulated time histories. The contours 
are construeted using the numerical values of the basis vectors plotted along the y-z plane for the 
heat shield section. Note that these charaeteristie patterns are similar to wave front propagation 
pulses. In addition to the patterns, numerieal values for the eontribution of the first 4 basis 
vectors to the overall response are shown. These values indieate that the basis eontributes 
13.5% to the overall response whereas basis number 2 eontributes 8.35%. Beeause the sum total 
eontribution from 4 bases is less that 33%, it should be apparent that many more basis veetors 
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participated in the overall response. As already mentioned, studying these patterns is very 
helpful to understand qualitatively the complexity of a pressure pulse, unfortunately, what is 
seen with a reduced set of sensors is far less. For example, figure 1 1 shows results for the first 
two basis veetors as seen through the 33 pressure sensors used in the simulated test. Certainly, it 
is extremely difficult to discern patterns but more importantly, the computed contribution is 
significantly different when compared to the full set of possible sensor loeations. 

7,3 Optimum Pressure Sensor Placement 

This section discusses an optimization based approach to place pressure sensors for this class of 
problems. The reader is reminded that this is strictly for illustration purposes because the 
approach was not implemented for testing. Nonetheless, the approach that follows recognizes 
that sensor placement must serve a dual role: 1) to calibrate the analytical model at locations 
where measurements exist, and 2) to allow one to make predictions at other locations. Often 
times these two roles are thought to be the same, but in fact that is not necessarily the case. For 
example, if there is local behavior captured in the model but not measured, then it is difficult to 
recover such behavior from measured data. In Section 5.2 an approach was discussed to perform 
response predictions using a reduced set of observation points and the analytical model. 
Effectively, the analytical model is used to provide interpolation functions to make predictions at 
other locations. Also discussed was the prediction error metric M3 to quantify the error incurred 
when conducting pressure predictions at other locations using only data from a reduced set of 
observation points. If sensors are placed properly, then predictions made throughout the surface 
will have minimum error. 

For the MLAS CM boiler plate test pressure sensors were placed (as shown in figure 2b) without 
regard to the ability to make predictions at other locations. Although placement using a heuristic 
approaches like the one used with MLAS can provide quality pressure data (as will be shown 
later), it hinders the ability to make predictions at other locations. To demonstrate this problem, 
a combinatorial optimization algorithm was used with simulated data to place 33 sensors. In here 
the sensor count was kept as in the test to maintain similar coverage. The optimization problem 
was tasked to choose 33 pressure sensor locations out of the 2148 possible locations while 
minimizing the prediction error metric M3. Figure 12a shows again the heat shield projected 
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onto y-z plane with the baseline sensors locations (o) and figure 12b shows the sensors placed 
using optimization (o). These results by themselves are not very revealing until one examines 
the prediction error. Figure 13a shows the normalized prediction error as a function of location 
when using responses from simulation at the baseline locations (figure 12a) to predict responses 
at all 2148 points. Certainly, at the 33 baseline locations the normalized prediction error is zero. 
Two aspects of the results should be apparent, 1) the normalized prediction error is very large 
(» 1) in some regions, and 2) there are clusters of regions where the error is small. Areas where 
the error is small are near measured points. However, in general the normalized error, defined as 
e = max(j^ - y) / max(j^) where je is the prediction and y is the exact value, is very large for many 
regions over the heat shield. In contrast. Figure 13b shows the normalized prediction error with 
sensors placed using optimization. For this case the normalized prediction error is three orders 
of magnitude smaller than results obtained using the baseline sensor set. Finally, to demonstrate 
the applicability of Eq. (1.3), figure 14 shows estimated pressures at two arbitrarily selected 
locations on the heat shield. Figure 14a shows time traces comparing pressure predictions 
(dashed) and the exact values (solid) at two different locations. Figure 14b show stars to indicate 
the location of these arbitrary points with respect to locations selected by the optimizer. Clearly, 
results shows that one is able to predict pressure responses at locations away from measured 
points. Nonetheless, the process can introduce significant numerical errors as evident in the 
prediction for pressure ID 5012903 in Fig. 14a. 


8 Concluding Remarks 

The paper discussed the application of newly developed calibration metrics to assess model 
adequacy in terms of probability of predicting the measured behavior and a multi-dimensional 
metric to measure model closeness. As a byproduct of implementing the metrics for model 
evaluation, a sensor placement approach was discussed and demonstrated to improve pressure 
predictions at points where measured data were not available. Unfortunately, the placement 
approach was developed after tests were completed and therefore could not be verified through 
test. Future testing should incorporate pre-test sensor placement as described in this paper to 
improve our ability to make pressure predictions at other locations. 
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To study the effects of variations in impact velocities and impact angles analytically, variations 
in these conditions were used to create a family of pressure pulse predictions. These results 
provided envelopes for the expected pressure predictions from the analysis as well as probability 
values. Data from two MLAS CM boilerplate tests were superimposed onto the analytical 
envelopes and it was determined that test results could not be enveloped by the analysis. Hence, 
test results cannot be reconciled with the analysis using the current model. Without additional 
data no improvements to the model are possible. 
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Table 1 Uncertainty parameters for T300-30 study 


No. 

Parameters 

Nominal 

Upper 

Lower 



Values 

Bound 

Bound 

1 

Pitch Angle (deg) 

30 

31 

29 

2 

Yaw Angle (deg) 

0 

0.4 

-0.4 

3 

Vertical Velocity (in/sec) 

300 

302 

293 


Table 2 Uncertainty parameters for T300-21 study 


No. 

Parameters 

Nominal 

Upper 

Lower 



Values 

Bound 

Bound 

1 

Pitch Angle (deg) 

21.5 

22.5 

20.5 

2 

Yaw Angle (deg) 

0 

0.4 

-0.4 

3 

Vertical Velocity (in/sec) 

300 

302 

293 
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Figure 1. MLAS test artiele and water impaet test at the U.S. Army Aberdeen test 



Figure 2b. Pressure sensor loeations 



Figure 3 a. MLAS FEM 



Figure 3b. Egg crate and local coordinate system 



Figure - Fluid Structure Interaction (FSI) Coupling Curve 


Eigure 4. Eluid Structure Interaction (ESI) coupling curve 



Figure 5. Pressure data uncertainty bounds for T300-30 (solid-blue) using 7 tests and 
uncertainty bounds for analysis (dashed-red) using 28 cases 
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Figure 6. Pressure data uncertainty bounds for T300-21 (solid-blue) using 7 tests and 
uncertainty bounds for analysis (dashed-red) using 35 cases 
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Figure 7. Orthogonality of measured and predicted impact shapes using pressure data for 


T300-30 


a) T300-21 Test8 
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Figure 8. Orthogonality of measured and predicted impact shapes using pressure data for 


T3 00-21 test case 8 at the top and test case 4 at the bottom. 



Figure 9. Augmented sensor set (+) and test sensor (o) loeations 
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Figure 10. Contour plots of the first 4 basis veetors using pressure data and 2148 


Heat 

Shield 

z-pos. 

(in) 


Basis No. 1 

Contribution 18.5997(%) 


Basis No. 2 
Contribution 11.6129(%) 



-0.0512 

120 



-0.0834 

110 




100 


. 

-0.116 

90 




80 



-0.148 


• 

• . 


70 


■ . • ♦ 



► 


-0.18 

60 




50 


' . * * 

-0.212 

40 



-0.245 

30 




on 




40 20 0 -20 -40 


0.3 

0.252 

0.205 

0.157 

0.11 

0.0B2 

0.0145 

- 0.033 

• 0.0806 

- 0.128 

- 0.176 

- 0.223 


Figure 11. Contour for first two pressure basis vector using 33 measured locations 
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Figure 12. Comparison of baseline pressure sensor placement and optimized placement 
using prediction error 



Figure 13. Comparison of normalized prediction error for pressure on the heat shield 
using baseline sensor placement and optimized sensor placement 
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Figure 14. Pressure predictions at unmeasured locations 


